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Abstract. Within the Green's function and equations of motion formalism it is possible to exactly solve a 
large class of models useful for the study of strongly correlated systems. Here, we present the exact solu- 
tion of the one-dimensional extended Hubbard model with on-site U and first nearest neighbor repulsive V 
interactions in the presence of an external magnetic field h, in the narrow band limit. At zero temperature 
our results establish the existence of four phases in the three-dimensional space (U, n,h) - n is the filling - 
with relative phase transitions, as well as different types of charge ordering. The magnetic field may dra- 
matically affect the behavior of thermodynamic quantities, inducing, for instance, magnetization plateaus 
in the magnetization curves, and a change from a single to a double-peak structure in the specific heat. 
According to the value of the particle density, we find one or two critical fields, marking the beginning of 
full or partial polarization. A detailed study of several thermodynamic quantities is also presented at finite 
temperature. 
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1 Introduction 

In physics exact solutions are of great importance since, 
in some cases, an approximation may introduce dramati- 
cally dominating errors, resulting in an incorrect descrip- 
tion of the phenomenon under study. Recently, within the 
Green's function and equations of motion formalism, we 
have shown that it is possible to exactly solve a large class 
of models useful for the study of strongly correlated sys- 
tems [HE]- By exactly solvable we mean that it is always 
possible to find a set of eigenenergies and eigenoperators 
of the Hamiltonian closing the hierarchy of the equations 
of motion. Thus, one can obtain exact expressions for the 
relevant Green's and correlation functions in terms of a 
finite set of parameters [3]. 

The aim of the present paper is twofold. First, we 
would like to further develop our previous work on the 
exact solution of the one-dimensional extended Hubbard 
model in the atomic limit (AL-EHM) [2], by extending it 
to a more general situation in which a finite magnetization 
may be induced by an external magnetic field. Secondly, 
the AL-EHM exhibits interesting features at low temper- 
atures, and one of the most fascinating characteristic fea- 
ture is that it shows magnetic plateaus. Interestingly, mag- 
netization plateaus have been predicted also for spin-one 
bosons in optical lattices [4]. By atomic limit, according 
to the conventional definition used in the literature, we 
mean the classical limit of the model, i.e., we set from the 
very beginning the hopping matrix elements Uj = 0. One 



may obtain different results by treating the problem with 
nonzero hopping and then taking the limit to zero. 

A uniform magnetic field, introduced through a Zee- 
man term, has dramatic effects both on the T = phase 
diagram and on the behavior of several thermodynamic 
quantities. We study the properties of the system as func- 
tions of the external parameters n, T/V, U/V and h/V 
(throughout the paper we set V — 1 as the unit of energy) , 
allowing for the on-site interaction U to be both repulsive 
and attractive. In fact, the parameter U can represent 
the effective interaction coupling taking into account also 
other interactions. Owing to the particle-hole symmetry, it 
is sufficient to explore the interval [0, 1] for the parameter 
n. The chemical potential is self-consistently determined 
as a function of the external parameters. We address the 
problem of determining the zero-temperature phase dia- 
gram in the (U, n, h) space and we find various phase 
transitions (PTs), as well as magnetization plateaus. At 
T = 0, for attractive on-site interactions, the magnetic 
field does not play any role if its intensity is h < \U\/2: 
the ground state is a collection of doublons (sites with 
two electrons of opposite spin), and there are no neighbor 
sites occupied. The magnetic energy is not strong enough 
to break the doublons. On the other hand, for strong re- 
pulsive on-site interactions, it is sufficient a small nonzero 
value of the magnetic field to have a finite magnetization. 
In the intermediate regions, the competition among U, V 
and h determines the phase structure. 
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Relevant thermodynamic quantities, such as the dou- 
ble occupancy, the charge and spin susceptibilities, the 
specific heat and the entropy are systematically computed 
both at T — and as functions of the temperature. For 
all values of the particle density, one finds a critical value 
(dramatically depending on U and n) of the magnetic field 
h a above which the ground state is ferromagnetic. In this 
state, every occupied site contains one and only one elec- 
tron, aligned along the direction of h. For strong repulsive 
on-site interactions, one finds h s = 0, i.e., the spin are po- 
larized as soon as the magnetic field is turned on. Further- 
more, for attractive on-site interactions and 0.5 < n < 1, 
one observes the existence of two critical fields, namely: h c , 
up to which no magnetization is observed, and h s , mark- 
ing the beginning of full polarization. This is analogous to 
the finite field behavior of the S — 1 Haldane chain [5J . 

The addition of a homogeneous magnetic field does not 
dramatically modify the framework of calculation given in 
Ref. [5] , provided one takes into account the breakdown of 
the spin rotational invariance. For the sake of comprehen- 
siveness, in Section [21 we briefly report the analysis lead- 
ing to the algebra closure and to analytical expressions 
of the retarded Green's functions (GFs) and correlation 
functions (CFs). The GFs and CFs depend on a set of in- 
ternal parameters, which can be determined by means of 
algebra constraints [UH], allowing us to provide an exact 
and complete solution of the one-dimensional AL-EHM 
in the presence of a magnetic field. In Section [3] we ana- 
lyze the properties of the system at zero temperature. We 
characterize the different phases emerging in the phase 
diagram drawn in the (U, n, h) space, by studying the be- 
havior of the chemical potential and of various local prop- 
erties (double occupancy, short-range correlation function, 
magnetization); Section [4] is devoted to the study of the 
finite temperature properties. Further to the study of the 
quantities analyzed in Section [3l we also present results 
for the charge and spin susceptibilities, the specific heat 
and the entropy. Finally, Sec. [5] is devoted to our conclu- 
sions and final remarks, while the Appendix reports some 
relevant computational details. 



2 The model 

The one-dimensional extended Hubbard model in the pres- 
ence of an external homogeneous magnetic field is de- 
scribed by the following Hamiltonian 

H = ~ 5ijli]c\i)c{j) + U ^n T (i)ri;(i) 



(1) 



n(i) = c^(i)c(i) is the charge density operator, c(i) (c^(i)) 
is the electron annihilation (creation) operator - in the 
spinor notation - satisfying canonical anti-commutation 
relations. 113(1) is the third component of the spin density 
operator, also called the electronic Zeeman term, 

= n]{i) - n^i) = c\{i) C] (i) - c\(i)c L (i). (2) 



Here we do not consider the orbital interaction with the 
magnetic field and we use the Heisenberg picture: i = 
(i, t), where i stands for the lattice vector R^. The Bravais 
lattice is a linear chain of N sites with lattice constant 
a. [i is the chemical potential; h is the intensity of the 
external magnetic field and U and V are the strengths 
of the local and intersite interactions, respectively. In the 
atomic limit, if one considers only first neighboring sites 
interactions, the Hamiltonian {!]) takes the form 



-fm(i) + UD(i) + -Vn(i)n a - hn 3 (i) 



(3) 



where D(i) is the double occupancy operator, defined as 
D(i) — n^(i)ni(i) — n(i)[n(i) — l]/2. Hereafter, for a 
generic operator <£(i), we define <£ a (i,t) = Ylj a ij^U^)t 
where onj is the projection operator over first nearest 
neighboring sites. 



2.1 Composite fields and equations of motion 

By taking time derivatives of increasing order of the fer- 
mionic field c(i), the dynamics generates other field opera- 
tors of higher complexity (composite operators). However, 
the number of composite operators is finite because of the 
recurrence relation 



[n a (i)] fc =X;^V(i)] r 



(4) 



which allows one to write the higher-power expressions 
of the operator n a (i) in terms of the first four powers. 
The coefficients Am are rational numbers satisfying the 

relations J2tn=i A ^ = 1 and A ^ = <W ( k = !> • ■ ■ ' 4 )' 
explicitly given in Ref. [B]. As a result, a complete set of 
eigenoperators of the Hamiltonian ([3]) can be found. 

Upon introducing the Hubbard operators = [n(i) — 
l]c(i) and rj(i) = n(i)c(i), one may define the composite 
field operator 



(5) 



where 



W(i)[n°(*)]V 



, ^ ) (i) = 



tfc,(i)[n«(i)] 
W(i)[n a (*)]V 



(6) 

where a = {T, I}- With respect to the case of zero mag- 
netic field [2] , the degrees of freedom have doubled, since 
one has to the take into account the two noncquivalcnt 
directions of the spin. By exploiting the algebraic prop- 
erties of the operators n(i) and D(i), it is easy to show 
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that the fields ip^ (i) and («) are eigenoperators of the 
Hamiltonian ([3|) [6]: 



at 

;W and are the 10 x 10 energy matrices: 



(7) 



e (0 = (A* o 



e 



(I) 



'e^ 



e 



(»/) 



(8) 



where 



r(fl 



/— /i-o-h 2V 

-fi-ah 2V 

-/M—crh 2V 

— /x— cr/i 

-3V #V — #V 







2V 

-H-ah+10vJ 



(9) 



,(>7) 



V 



M- 


ct/i 2V 








\ 





U—fj>—<rh 


2 V 














U — ti — crh 


2V 














U — [i — crh 


2\/ 





-3V 


fv 


-fv 


[/-p-crh+lOV/ 



The eigenvalues of the matrices £^ and are 



CS) N 



(€) 



E 



In) 



where 



E® =-»-ah+(p-l)V, 
E$ = U-n-<Th+(p-l)V, 



(10) 



(11) 



(12) 



with p = 1, . . . , 5. Although at the level of equations of 
motion the two fields tp^(i) and tfj^ii) are decoupled, 
they are indeed coupled by means of the self-consistent 
equations which determine the correlators appearing in 
the normalization matrix, as shown in Subsection 12.31 



2.2 Retarded Green's and correlation functions 

The knowledge of a complete set of eigenoperators and 
eigenenergies of the Hamiltonian ([3|) allows one to exactly 
determine the retarded thermal Green's function 

G<°>(i - t') = 9(t - OW (o) (i.*)^ ot (i,*')}>- (13) 

In the above equations, a = £, r\ and (• ■ • ) denotes the 
quantum-statistical average over the grand canonical en- 
semble. By means of the field equations (JTJ) , it is easy to 
show that the Green's function satisfies the equation 



[u-£t , )]GW(w)=/W 1 



(14) 



where I^ a > is the normalization matrix 

/ (o M{V> (o HM),V> (o)t (M)}>, 



(15) 



whose expression in discussed in the Appendix. The solu- 
tion of Eq. (HU) is 
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p=i w - E i 

(a,p) 



(a) 



iS 



(16) 



where the spectral functions <Jjj,v can be computed by 
means of the formula 



n {a,p) _ o(a) (Q) V a) 



(17) 



J?( a ' is the 10 x 10 matrix whose columns are the eigenvec- 
tors of the matrix gW. Explicit expressions of the matrices 
= QW) are reported in the Appendix. The correlation 
function 

C {a) {t-t') = (ip^(i,t)ip {a) \i,t')) 

1 f +oo (is) 



(27T) J- 

can be immediately computed from Eq. by means of 
the spectral theorem, and one finds 
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CW(a,)=7r£ 
p=i 



1 - lauJ] | ~y 



a MS[u-E^], (19) 
where (3 = 1/ksT. At equal time, the CF is given by 

10 

C (Q) = (^ (Q) (i)V (a)t W) = ^53 T i 0)<j(0 ' P) ' ( 2 °) 



P =i 



where Tp°' = l + t&nh(E p a ^ /2ksT). Upon introducing the 
parameters 



n = (n(i)), 

m = (n 3 (i)), - (21) 

C W = ([n°(i)n, ^ (P) = d«3(i)KWF), 



A (p) = 5 (n(i)[n«(i)P>, 



the elements of the normalization matrix can be written 
as 

I® = K<r-V-\<r- 1 '> + n<r- 1 \ 
I® + z = k<*-V-\<*- 1 '> 

tW _ A (p-i) , (p-i) 



(22) 



where and p = 1, . . . , 5. All the other matrix elements can 
be easily computed by means of the recurrence relation 
©■ 
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2.3 Self-consistent equations 

The previous analysis shows that the complete solution of 
the model requires the knowledge of the following 13 pa- 
rameters: fi, m, , . . . , k^ 4 ) , , . . . , A^ 4 ) , 
These quantities may be computed by using algebra con- 
straints and symmetry requirements [TJ[3] • By recalling the 
projection nature of the Hubbard operators and T](i), 
it is straightforward to verify that the following algebraic 
properties hold 



e t ( i )n(z)=0, 
^(i)D(i) = 0, 



rf{i)D{i) = 0. 



(23) 



These are fundamental relations and constitute the ba- 
sis to construct a self-consistent procedure to compute 
the various parameters of the model. Upon splitting the 
Hamiltonian ([3]) as the sum of two terms: 



H — H + Hi, 
Hi = 2Vn[i)n a {i), 



(24) 



the statistical average of any operator O can be expressed 
as 



(0) = 



(O 



-mi 



>0 



-I3H, 



>0 



(25) 



where (• • • )o stands for the trace with respect to the re- 



duced Hamiltonian Hi 



o- 



Tr{--- 



-pH 



}/Tr{e-^ } 



The Hamiltonian Hq describes a system where the original 
lattice has been reduced to a central site and to two un- 
connected sublattices. Thus, in the .ffo-representation, the 
correlation functions connecting sites belonging to discon- 
nected sublattices can be decoupled. Within this scheme, 
the unknown parameters can be written as functions of 
only two parameters X\ — (n a (i)) and X 2 — (D Q (i))o, 
in terms of which one may find a solution of the model. 
By exploiting the translational invariance along the chain, 
one can impose (n(i)) = (n a (i)) and (D(i)) = (D a (i)), 
finding, thus, two equations allowing one to determine X\ 
and X2 as functions of \x: 



F 1 (X 1 ,X 2 ; t i) = 0, 
F 2 (X 1 ,X 2]fJ ,) = 0. 



(26) 



The chemical potential /i can be determined by means of 
the equation 

n = F 3 (X 1 ,X 2 ;n). (27) 

The explicit expressions of the functions F\ , F 2 , and F 3 are 
given in the Appendix. Equations ([26]) and (|27p constitute 
a system of coupled equations allowing us to ascertain the 
three parameters [i, X\ and X 2 in terms of the external 
parameters of the model n, h, U, V, and T. Once these 
quantities are known, all the properties of the model can 
be computed. 



3 The zero-temperature phase diagram 

In this Section we derive the phase diagram of the AL- 
EHM in the (U, n, h) space. By numerically solving the set 




Fig. 1. The phase diagram in the space (U, n, h) at T — and 
V = 1. The four different phases are characterized by differ- 
ent distributions of the particles. Phase (a) is characterized by 
zero magnetization and by only non-neighboring doubly occu- 
pied sites. In the phase (b) the electrons singly occupy non- 
neighboring sites with all the spins parallel to h. In the phase 
(c) one finds both singly polarized and doubly occupied non- 
neighboring sites, whereas in the phase (d) neighboring sites 
can be occupied but no doublons are to be found. 



of equations (|26|) and (|27j) , we study the T = behavior 
of relevant physical quantities. This investigation allows 
us to envisage the distribution of the particles along the 
chain, for different densities, as well as the magnetic prop- 
erties of the system. The results obtained are displayed in 
Fig. [1] one recognizes the already known [3] four different 
phases found at h = extending in the h direction. 

The phase structure is determined by the three com- 
peting terms of the Hamiltonian: the repulsive intersite 
potential (disfavoring the occupation of neighboring sites) , 
the magnetic field (aligning the spins along its direction, 
disfavoring thus double occupancy) and the on-site poten- 
tial, which can be either attractive or repulsive. According 
to the values of these competing terms, one may distin- 
guish the different phases, characterized by different val- 
ues of the double occupancy, the chemical potential, and 
of the parameters defined in Eq. (|2~lj) . 

Phase (a) is observed in the region < n < 1 for at- 
tractive on-site potential (U < 0) and h < \U\/2; it is 
just a continuation along the h axis of the phase observed 
at h = f2] and is characterized by a zero magnetiza- 
tion m = 0. The parameters Xi and X 2 take the values 
X\ = 2X 2 = 2n/(2 — n). The chemical potential takes 
the value U/2 for n < 1, whereas at n = 1 (half fill- 
ing) it jumps to the value /1 = 2V + U/2, as required 
by particle-hole symmetry. The distribution of the parti- 
cles in the phase (a) is shown in Fig. [21 where we report 
just one possible configuration. The attractive on-site po- 
tential favors the formation of doubly occupied sites. At 
the same time, the nearest-neighbor repulsion V disfa- 
vors the occupation of neighboring sites. This distribution 
of the electrons is confirmed by the expectation values 
D = n/2 and = 0. When n< 1, there is no ordered 
pattern in the distribution of the particles, whereas, for 
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n=0.25 



n=0.25 



n=0.5 

n=0.75 

n=l 



n=0.5 

n=0.75 

n=l 



Fig. 2. Distribution of the particles along the chain by varying 
the particle density at T = 0, V = 1, U < and h < \U\/2. 
White and black squares denote empty and doubly-occupied 
sites, respectively. 



n = 1, one observes the well-known checkerboard distri- 
bution of doubly-occupied sites [HE] . The phase (6) is ob- 
served only for particle densities equal or less than quarter 
filling (n = 0.5), in the regions (U > 0, V7i) and (U < 0, 
h > \U\/2). In the latter region, there is a competition 
between the attractive on-site potential which favors the 
creation of doublons and the magnetic field which favors 
the alignment of the electrons. When h > \U\/2, it is en- 
ergetically convenient to singly occupy non-neighboring 
sites. Similarly, for U > 0, the repulsion between elec- 
trons on the same site and on neighboring sites (V > 0), 
leads to a scenario where the double occupancy, as well as 
the short-range correlation functions \( k ' , vanishes in the 
limit T — > 0. One observes that the electrons tend to singly 
occupy non-neighboring sites with all the spins parallel to 
h, leading to a finite magnetization m = n. For particle 
densities less than quarter filling, there is a cost in energy 
to add one electron which is proportional to the intensity 
of the magnetic field: /z = —h. At n = 0.5, the chemical 
potential jumps to the value \i = V — h and a long-range 
order is established: one observes a checkerboard distribu- 
tion of singly occupied polarized sites, as evidenced in Fig. 
[3] In this region one finds X\ = n/(l — n) and X 2 = 0. 

For particle densities greater than quarter filling, the 
phase diagram is richer in the plane (U,h): one observes 
three phases by varying U (h) and keeping h (U) fixed, as 
evidenced in Fig. [TJ Phase (a) has been discussed above. 
The phase (c) is observed for 0.5 < n < 1 and in the 
regions (0 < U < 2V, h < V - U/2) and (U < 0, 
\U\/2 < h < V + \U\/2). The chemical potential takes 
the value U + h for 0.5 < n < 1 and 2V + U/2 at n = 1. In 
this region the intersite interaction dominates: the mini- 
mization of the energy requires the extra electrons above 
n = 0.5 to occupy non-empty sites. Thus, in the thermo- 
dynamic limit, there can be singly polarized and doubly 
occupied sites but no neighboring sites: (n(i)[n a (i)] k ) = 0. 
By increasing n, the number of doubly-occupied sites in- 
creases linearly, D = n — 1/2, while the magnetization 
decreases as m = 1 — n. The distribution of the parti- 
cles in the region (c) is drawn in Fig. [3] where we rep- 
resent just one possible configuration. The ground states 
of the Hamiltonian are checkerboard configurations where 
empty sites alternate with sites occupied by either two 
particles or one particle with polarized spin. At half fill- 
ing, one observes a checkerboard distribution of only dou- 
bly occupied sites [7j. Also the phase id) is observed for 
both attractive and repulsive on-site interactions. In par- 
ticular, when U > 2V the on-site interaction dominates 
over the nearest-neighbor repulsion V. The minimization 



Fig. 3. Distribution of the particles along the chain by varying 
the particle density at T — and V = 1. For < n < 0.5, 
this distribution is observed both for {U < 0, h > |J7|/2} 
and {U > 0, V/i}. For 0.5 < n < 1, it corresponds to the 
regions {U < 0, \U\/2 < h < V + \U\/2] and {0 < U < 2V, 
h < V — U/2}. White and black squares denote empty and 
doubly-occupied sites, respectively. The squares with the arrow 
denote singly occupied sites with field-aligned spins. 



of the energy requires the electrons not to be paired and 
allows for the occupation of neighboring sites: D = and 
A^ 1 ) =rt — l/2. The combined action of U and h predomi- 
nates over V, leading to the same distribution in the region 
(0 < U < 2V, h > V-U/2). For U < 0, a strong magnetic 
field (h > V+\U\/2) dominates both the on-site and inter- 
site potentials, resulting on the absence of doublons and 
the possibility to have neighboring sites occupied. In this 
phase, one finds X\ = l and X 2 = 0. For n < 1 the energy 
necessary to add one electron is /x = 2V — h. At n = 1 the 
chemical potential jumps to the value \i — 2V+U/2. Since 
in this state every occupied site contains one and only one 
polarized electron, the ground state is ferromagnetic and 
the magnetization is m = n, as evidenced in Fig. [4j 



MtM,MtM,MtM,|t|tM,|t|tM,|t|t|,MtM,Ht|t|,MTm n=l 

Fig. 4. Distribution of the particles along the chain by varying 
the particle density at T = and V = 1. This distribution is 
observed in two different parameters regions, namely {U < 
2V, h > V + U/2} and {U > 2V,V/i}. White squares denote 
empty sites, whereas the squares with the arrow denote singly 
occupied sites with field-aligned spins. 



These peculiar distributions of the electrons along the 
chain by varying the magnetic field give rise to the forma- 
tion of plateaus in the magnetization curves. By increas- 
ing the magnetic field, there are plateaus whose starting 
points depend on the particle density, as well as on the 
on-site potential: one identifies two critical values of the 
magnetic field. The nonzero magnetization can either be- 
gin from h — or from a finite field. h c denotes the start- 
ing point of a nonzero magnetization, whereas h s denotes 
the value of the magnetic field when it reaches saturation. 
The results for the magnetization m(h) are shown in Fig. 
[5] The values of h c and h s can be inferred from the previ- 
ous analysis. Making reference to Fig. [5] for the different 
regions of n and U, one has: h c — h s — \U\/2 (Fig. [5^), 
h c = \U\/2 and h s = V + \U\/2 (Fig. [5J>), h c = and 
h s = V- U/2 (Fig. [5fc), and h c = h s = (Fig. EH). 
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(d) 



V=1 
U>0 
0<n<0.5 



f -1 
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T=0.1 

- 7=0.05 
7=0 



T=0.2 
■7=0.1 

- 7=0 05 

- 7=0 



Fig. 5. The magnetization m as a function of the magnetic 
field h at V = 1 and decreasing temperatures. 



The so-called metamagnetic behavior is clearly seen: 
at low temperatures the magnetization begins to show a 
typical 5-shape which becomes more pronounced by fur- 
ther decreasing the temperature. At T = one, two or 
three plateaus (to = 0, to = 1 — n and m — n) are ob- 
served, according to the values of the external parame- 
ters. These results are similar to the ones obtained in a 
one-dimensional spin-1 antiferromagnetic Ising chain with 
single-ion anisotropy [UGH- 

In Figs. [6] we plot the double occupancy and the chem- 
ical potential as functions of h for U = —1, n = 0.25, 
0.75 and different temperatures. As the temperature de- 
creases, D exhibits sharp jumps, whereas fi shows discon- 
tinuity (ies). In particular: (i) for all values of the particle 
density, fi takes the constant value fi = U/2 in the re- 
gion U < and h < \U\/2, corresponding to h s (h c ) for 
n = 0.25 (0.75); (ii) for < n < 0.5 and h > h s , fi 
decreases with the law \i — —h. On the other hand, for 
0.5 < n < 1, fi first increases with the law fi = U + h 
until h s = V + \U\/2. Further increasing h one observes 
a decrease of the chemical potential following the law 
fi = 2V — h. As it is evident from Figs. El at T = and for 
repulsive on-site interactions, the double occupancy shows 
a one or a two step-like behavior according to the parti- 
cle density. The steps occur at the values of the magnetic 
field where PTs are observed. In closing this Section, it is 
worthwhile to mention that at T = and h = all the 
phases (with the exception of the case n = 1 and U < 2V) 
exhibit a macroscopic degeneracy growing exponentially 
with the volume of the lattice. A nonzero magnetic field 
can lift the degeneracy of the ground states at quarter and 
half filling, even for repulsive on-site interactions. 



Fig. 6. The chemical potential and the double occupancy at 
V = 1, U = — 1 and different values of the temperature, for 
n = 0.25 and n = 0.75. 



system as a function of the parameters n, T, U and h, and 
again we take V = 1 as the unit of energy and we set the 
Boltzmann's constant kg = 1- 



4.1 Thermal properties 

In Fig. [7] we show the behavior of different quantities as 
a function of T for n = 0.5, U = —1 and various values 
of h in the vicinity of the critical point h s = \U\/2 = 0.5. 
These behaviors are not peculiar to the given values of 
the parameters, but are always observed in the neighbor- 




4 Finite temperature 

In this Section we shall investigate the finite temperature 
properties of the AL-EHM. We study the behavior of the 



Fig. 7. (a) The double occupancy, (b) the chemical potential, 
(c) the magnetization m, and (d) the charge susceptibility as 
functions of the temperature for V = 1, U = —1 and n = 0.5 in 
the neighborhood of h s — 0.5, where the transition from phase 
(a) to phase (b) occurs. 
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Fig. 8. (a) The spin susceptibility as a function of the tem- 
perature for V — 1, U = —1, n = 0.75 in the neighborhood of 
h = \U\/2. (b) The entropy S as a function of the temperature 
at V = 1, U = 1, n = 0.5 and different values of h. 



hood of all critical values of h when the system goes from 
one phase to another. In the high temperature regime the 
double occupancy, the chemical potential and the mag- 
netization decrease by increasing T, whereas for T — > 
they tend to two different constants, as shown in Figs. [7l 
At h = h s and T — there is a phase transition: for the 
considered values of the external parameters, one finds 
that the double occupancy jumps from zero to n/2, the 
magnetization from zero to n and the chemical potential 
from U/2 to V — h. Also the charge susceptibility, de- 
fined as \c — N n 2 + ks T dn/d^i, behaves differently at 
low temperatures according to the value of h. For h < h s 
the charge susceptibility increases by decreasing T and 
tends to Xc = n(l — n){2 — n) at T = 0. For h > h s the 
charge susceptibility decreases by decreasing T and tends 
to \c = n(l — n)(l — 2n) at T = 0. In the particular case 
shown in Fig. [TJi, Xc vanishes for T — > when h > h s : the 
system is in a charge ordered state [2] ■ The non- vanishing 
of the charge susceptibility in the limit T — > corresponds 
to a non-ordered ground state. For h = h s and T = there 
is a phase transition and the charge susceptibility exhibits 
the discontinuity Ax c — n{l — n 2 ). 

Another signature of the T — transition in the neigh- 
borhood of the critical values of the magnetic field is pro- 
vided by the behavior of the spin susceptibility, defined 
as Xs — dm/dh. As an example, in Fig. we plot Xs 
for U = — 1, n — 0.75 and different values of h around 
h c = \U\/2. In the high temperature regime the spin sus- 
ceptibility decreases by increasing T. In the low tempera- 
ture regimes, the spin susceptibility exhibits a peak at a 
temperature T±, then decreases going to zero at T = 0. 
By approaching the value h = h c (both from below and 
above) the position T\ of the peak moves towards lower 
temperatures; at the same time the height h\ of the peak 
increases. At h — h c the spin susceptibility diverges at 

r = o. 

Recently, we have discuss how the entropy may be 
computed when the ground state is degenerate [2] . In this 
case, the zero-temperature entropy S(0) is non-zero and, 
if the system is not confined in one of the possible phases, 
is not even constant but depends on the external param- 
eters. The magnetic field can remove the h = ground 
state degeneracy, responsible of a finite zero-temperature 



Fig. 9. The specific heat C as a function of the temperature 
for V = 1, U = —1, n — 0.5 and for different values of h. 

entropy. This can be clearly seen in Fig. [SJd where we plot 
the entropy as a function of the temperature at U = 1 
and n = 0.5 for increasing values of the magnetic field. 
As soon as the magnetic field is turned on, the entropy 
vanishes in the limit T — > 0. 

The study of the specific heat further enlightens the 
influence of the magnetic field on the thermodynamic be- 
havior of the system. The specific heat is given by C = 
dE/dT, where the internal energy E can be computed as 
the thermal average of the Hamiltonian ([3|) and it is given 
by E = UD + 2V\^ — hm. As an example of the charac- 
teristic behavior of the specific heat by varying h, in Fig. 
HO we plot C as a function of the temperature at V = 1, 
U = —1, and n = 0.5, for several values of the magnetic 
field. With the exception of the critical value h = h s , one 
observes an exponential activation of C with a pronounced 
peak, whose position depends on h. Furthermore, the low 
temperature peak observed at n < 1 for pertinent values 
of the on-site potential when h = 0, tends to disappear 
by approaching h s , both from above and below. The po- 
sition of the peak moves towards lower temperatures and 
vanishes exactly at h — h s . Of course, also the low temper- 
ature peak with exponentially increasing height (observed 
at half-filling in the limit U — > 2V at h — [2]) survives 
only for small values of h, but disappears at h = h s . On 
the other hand, for higher values of the magnetic field one 
always observes a double peak structure in the specific 
heat for all values of the particle density. In this region, 
the low temperature peak remains constant by varying h, 
whereas the high temperature broad peak moves towards 
higher temperatures by increasing h. 

Interestingly, a change from a single to a double-peak 
structure by varying the magnetic field has also been ob- 
served in one-dimensional ferrimagnets with alternating 
spins [TO] . 



4.2 Magnetic properties 

A useful representation of the phase diagram is obtained 
by plotting the thermodynamic quantities as functions of 
the magnetic field h. Indeed, this representation consti- 
tutes another way to detect the zero-temperature tran- 
sitions from thermodynamic data. In fact, at low tem- 
peratures, all the thermodynamic quantities that we have 
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Fig. 10. The charge susceptibility \ a as a function of the 
magnetic field at different temperatures for V = 1, n = 0.75 
and for (a) U = -1, (b) (7 = 1, and (c) U = 2.5. (d) Xc as a 
function of the particle density for V = 1, U = — 1, T = 0.01 
and different values of h. 



Fig. 11. The spin susceptibility as a function of the magnetic 
field for V = 1, T = 0.1 and (a): U = -1, (b): U = 1, (c): 
t/ = 2.5. 



investigated present minima, maxima, or discontinuities 
in the neighborhood of the critical values of h at which 
zero-temperature PTs occur. 

As a first example, in Figs.fTOk-c we plot Xc as a func- 
tion of the magnetic field, for n = 0.75 and for different 
values of the on-site potential (U = —1,1, 2.5). In the high 
temperature regime, the charge susceptibility is always a 
decreasing function of h. In the limit T — > 0, one observes 
two different behaviors, according to which phase the sys- 
tem is in: \ c tends to a constant when the system is in a 
non-ordered phase, whereas it presents minima in corre- 
spondence of the critical values of the magnetic field, as 
evidenced in FigfTUJ A similar behavior of the charge sus- 
ceptibility is observed also for other values of the particle 
density. In Fig. fTOt l we plot Xc as a function of the filling 
n at V — 1 and U = — 1 and for two values of h (h — 0.45 
and h — 0.55, below and above h cr u — 0.5, respectively). 
To simplify notation, here we shall use h cr i t to indicate ei- 
ther h c and h s . One immediately sees that, for T — > and 
h < h cr it-, the charge susceptibility increases by increasing 
n and has a maximum at quarter filling. Further increas- 
ing n, Xc decreases and vanishes at half-filling, where, at 
fixed on-site potential, a zero-temperature PT occurs by 
varying n, from a non-ordered to a checkerboard distri- 
butions of doublons (see Fig. [5]). When h > h crit , at low 
temperatures, Xc has a double peak structure with two 
maxima around n = 0.25 and n = 0.75 and two minima 
around n — 0.5 and n = 1. This implies the presence of a 
charge ordered state also at n = 0.5, as evidenced by the 
charge distribution shown in Fig. [4] 

In Figs. [TTk-c we plot the spin susceptibility as a func- 
tion of the magnetic field at T = 0.1 for three represen- 
tative values of U (U = — 1, [7=1 and U = 2.1, respec- 
tively) and for different values of the filling (n = 0.25, 0.5, 



0.75, 1). The spin susceptibility diverges in correspondence 
of the values h cr it at which one moves from one magneti- 
zation plateau to the other. For low values of the magnetic 
field and attractive on-site interactions - corresponding to 
Fig. II lb - the spin susceptibility tends to vanish at low 
temperatures for all values of the filling: at T = all elec- 
trons are paired and no alignment of the spin is possible. 
By increasing h, the magnetic excitations break some of 
the doublons inducing a finite magnetization: Xs has a 
peak, then decreases, the system having entered the suc- 
cessive magnetic plateau. If 0.5 < n < 1 then another 
peak is observed, corresponding to the second jump of the 
magnetization when h reaches the saturated value h s . On 
the other hand, for repulsive on-site interactions, a very 
small magnetic field induces a finite magnetization (with 
the exception of n = 1 when U < 2V): Xs has a max- 
imum at h = and decreases by augmenting h, unless 
another transition line is encountered, as it happens for 
0.5 < n < 1 and < U < 2V. 

In Figs. [12] we plot the entropy as a function of the 
magnetic field for relevant values of the particle density 
n. The entropy presents maxima in the neighborhood of 
the values of h at which one observes zero-temperature 
PTs, and 5(0) has discontinuities right at those values. 
On the other hand, for sufficiently strong magnetic fields, 
the entropy becomes rather insensitive to variations in h. 
When n < 0.5, the entropy presents a peak around h a 
which becomes more pronounced as the temperature de- 
creases. In the region 0.5 < n < 1, at low temperatures, 
one finds two peaks at h c and h s , respectively. Although 
at half filling there is a zero-temperature transition at h s , 
one always finds 5(0) = since the relative ground states 
are non-degenerate. 
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Fig. 12. The entropy 5 as a function of the magnetic field for 

V = 1, U = -1, T = 1, 0.5, 0.2 and (a): n = 0.25, (b): n = 0.5, 
(c): n = 0.75, and (d): n = 1. 

5 Concluding remarks 

We have evidenced how the use of the Green's function 
and equations of motion formalism leads to the exact so- 
lution of the one-dimensional AL-EHM limit in the pres- 
ence of an external magnetic field. We provided a system- 
atic analysis of the model for nearest-neighbor repulsion 

V by considering relevant thermodynamic quantities in 
the whole space of the parameters n, T, h and U (hav- 
ing chosen V — 1 as the unity of energy) . This study has 
shown that, at zero-temperature, the model exhibits phase 
transitions for specific values of n, h and U. In particular, 
we have identified four phases in the (U,n,h) space and 
PTs are observed at the borders of these phases. Various 
types of long-range charge ordered states have been ob- 
served: (i) at half-filling for U < 2V and h < V - U/2, 
a checkerboard distribution of doubly-occupied sites; (ii) 
at quarter-filling for U > 0, as well as for U < and 
h > \U\/2, a checkerboard distribution of polarized singly 
occupied sites; (iii) for 0.5 < n < 1 an ordered state 
with alternating empty and occupied sites in the regions 
(U < 0, \U\/2 < h < V + \U\/2) and (0 < U < 2V, 
h>V- U/2). 

We derived the phase diagram in the space (U, n, h) at 
T = by computing several quantities whose behaviors 
is also useful to characterize the distribution of the elec- 
trons on the sites of the chain. When plotted as functions 
of h, the chemical potential, the double occupancy as well 
as the magnetization, show discontinuities where PTs oc- 
cur. We identified the values of the critical fields h c and 
h Sl defining, according to the value of the particle density, 
the beginning point of nonzero magnetization and the sat- 
urated magnetization field, respectively. Furthermore, the 
presence of the magnetic field may dramatically modify 
the behavior of several thermodynamic quantities. For in- 
stance, the charge susceptibility tends to different values 
in the limit T — > 0, depending on how strong is the mag- 



netic field. At quarter and half filling, a strong magnetic 
field can also remove the macroscopic degeneracy of the 
ground states: as a result, the zero-temperature entropy 



A Computational details 

Firstly, we shall provide some more details on the matrices 
used in the computational framework adopted in Sec. [2] 
To begin with, we note that, due to the spin conservation, 

the normalization matrix 1^ = ({i/> (a) (i, i), i/> (a)t (i,t)}> 
is block diagonal: 



/TT,(o) o 
/U>(°) 



(28) 



The time translational invariance requires the m-matrix 
m = {{i§i^{i, t), ip^(i, t)}) to be symmetric. This require- 
ment leads to find that several elements of the normaliza- 
tion matrix are not independent, and one needs to com- 
pute only the matrix elements {h,p}, and {Ig. p +g}, where 
p = I, ... ,5. By recalling the basic commutators 



{£,(*),£(*)} = [l-«-a(i)]. 

{VtT(i),Tit(i)} = n^ a {i), 

one easily obtains Eq. (f22|) . 

The matrices Q^> and have the expression 



(29) 



(\ 2 4 1 (§) 4 

2 3 1 (|) 3 

2 2 1 (§) 2 

2 1 1 (§) 

2° 10 (§) 



(j) 
* 



1 2 4 1 

2 3 1 

2 2 1 

2 1 1 

^0 2° 1 



(I) 
(I) 1 



(30) 



As it has been shown in Sec. [H the Green's and the 
correlation functions depend on 13 parameters: /x, m, k^ 2 \ 
K W, A« A( 2 ), AW, AW, ttW, tt^, tt^, tt^. By 
recalling the expression of Hj [see Eq. (|24p] . and by using 
the recurrence rule (J3J, it is easy to show that 



H ' = 1 + J2 H^fm + D(i)g m ] K«p, (31) 



rn—1 



where 



oo 1 
n— 1 

A ^ Tjl 



(32) 
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The coefficients f m and g m are polynomials of K — e " v 
given by 



h = q(K - i) 2 (hk 2 - UK + 35) 

f 3 = ~(K-l) 3 (3K-5), 



(33) 



and 

9i 

92 
93 



6 



(X - l) 5 (3K 3 + 15K 2 + 29K + 25) 



-(K - l) 2 (UK 6 + 22K 5 - 23K 4 - 68K 3 
21K 2 + 138K - 35) , 

--{K - l) 3 (3K 5 + 9Jf 4 + 4K 3 
3 

12K 2 - 21K + 5) , 



(34) 



54 = r{K - l) 4 (if 4 + 4if 3 + 6if 2 +4K - 1) 
3 



Now, by taking the expectation value of Eq. ([31]) with 
respect to Hq, and exploiting to the properties of the Hq- 
representation, one finds 



To = {e- 0H % = 1 + (n(i))o E -UKWDo 

771—1 

4 

+ (j)(»))„^ Sm (K(i)P) , 



(35) 



In order to compute the quantities (n(i))o and (-D(i))o, 
one may use the equations of motion 



[£ a (i),H ] =-{ i i + ah)^{i), 
[r) a (i), H ] = -(fi + crh- U)rj,y(i), 



(36) 



leading, for a homogeneous phase, to the following expres- 
sions for the CFs: 



1 - (n_ <T (i)) n 
1 _|_ e P(n+vh) ' 



(37) 



Recalling that + r^J. = 1 - n CT and 77^77* = 7i_ CT - 
TijTij, from Eq. ([37]) one obtains 



0(fl+2h) + e P(2fi+h-U) 



B{ = (ni(i)) 



3 /3/i _|_ + e /3(p+2h) + e P(2^+h-U) ' 

e /3M _|_ e f)(2^+h-U) 
o/3h _|_ e /3 M _|_ e /3(/i+2h) _|_ e p(2/i+h-U) ' 



Si = bJ + si 



B 2 = <D(i)) 



; e /3M(i + e 2/3ft + 2 e^^+' 1 - c/ )) 



B 3 = B\ - b[ = 



j/3/i _)_ e /3/j _)_ e /3(ii+2h) _|_ e /3(2p+/i-[/) ' 



(38) 

Because of the properties of the 7?o- re P resen tation, by us- 
ing Eq. (|3T[) and by means of the relations 



n a (i)n(i) = n a (i) + D(i), n 3 (i)n(i) = n 3 (i), 
n a {i)D{i)=D{i), n 3 {i)D(i) = 0, 



(39) 



one finds 



T ° = 1 + J2 ( B *f™ + B2g m ){{n a (i)} m )o, (40) 



and 



1 4 

n = ^r{B 1 + J2 [( fi i + 2B 2)/ m + 2B 2 g m ]{[n a (i)] m ) }, 



To 



m— 1 
4 



{i+E/ m (K(7)Do}, 



m— 1 
4 



D = t^ 1 + E( 2 /™ + ffmXK«] m >o}. 

jo 



Upon defining the two parameters 

*i - (n°(i))o, 
X 2 = (D a (i)) , 



(41) 



(42) 



and by exploiting the translational invariance along the 
chain, (n(i)) — (n a (i)) and (D(i)) — (D a (i)), one obtains 
two equations allowing one to determine X\ and X 2 as 
functions of namely: 

4 

F 1 =X 1 -B 1 +J2 (Bifm + B 2 g m )([n a (t)] m+1 ) 

m=l 

4 

- E [(^i + 2S 2 )/ m + 2B 2 g m ]([n a (t)} m ) = 0, 

m=l 

4 

F 2 = X 2 - B 2 + (Bifm + B 2 g m )(D a (i)[n a (i)] m ) Q 

m=l 

4 

-i? 2 E(2/ m +ff m )(K«] m )o = 0. 



m— 1 



(43) 
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Upon denning S = (1 - X Y + X 2 ) + {X 1 - 2X 2 ) (1 + aX 1 + 
a 2 X2)+X2(l + dXi + d 2 X2), the chemical potential /x can 
be determined by means of the equation 



F 3 = - 2X 2 )(1 + aX x + a 2 X 2 ) 

+ 2X2(l + dX 1 +d 2 X 2 ) = n. 



(44) 



After some lengthy but straightforward calculations, one 
finds that Eqs. (|43|) can be rewritten as 



Xi = 2e /3AI cosh(/3/i)(l -X x - dX 2 )(l + aX x + a 2 X 2 ) 
+ e ^- u) [2 + (d- l)Xi - 2dX 2 ]{l + dXi + d 2 X 2 ), 

X 2 = e t3 ^~ u) [l + dXi - (2d + 1)X 2 ](1 + dXi + d 2 X 2 ) 
- 2e^ cosh(/3/i)(l + d)X 2 {l + aX x + a 2 X 2 ). 

(45) 

As a consequence, one also gets: 

D = — Jf 2 (l + dXi + d 2 X 2 ), 
1 



= — tanh(/3/i)(Xi - 2A 2 )(1 + aXi + a 2 X 2 ) 



= tanh(/3/i)(n - 2D), 
A« - ^^[(^1 - 2X 2 )(Xi + 2aX 2 ) 
+ 2KX 2 (X 1 + 2dX 2 )}, 



(46) 



where Q = + 2GK + K 2 ) 2 - 8GK(K - l) 2 . Upon, 
inserting Eq. (29|) in Eqs. (44|) and ((46]) one finds 



1, 



D = 



1 



A (1 



2[l + G(l-a 2 A 2 ) 2 ]' 
1 - 2a 2 AT 2 + (1 — AK ■ 



(50) 



2if 2 + K 4 )X 2 



2[l-2a 2 A 2 + a 2 (l + X 2 )Asr|] 



References 

1. F. Mancini, Europhys. Lett. 70, 484 (2005); Condens. Mat- 
ter Phys. 9, 393 (2006). 

2. F. Mancini, F. P. Mancini, Phys. Rev. E 77, 061120 (2008). 

3. F. Mancini, A. Avella, Adv. Phys. 53, 537 (2004). 

4. A. Imambekov, M. Lukin, E. Demler, Phys. Rev. Lett. 93, 
120405 (2004). 

5. F. D. M. Haldane, Phys. Lett. A 93, 464 (1993); Phys. Rev. 
Lett. 50, 1153 (1983). 

6. F. Mancini, Eur. Phys. J. B 47 527, (2005). 

7. R. A. Bari, Phys. Rev. B 3, 2662 (1971). 

8. X. Y. Chen, Q. Jiang, W. Z. Shen, C. G. Zhong, Jour. Magn. 
Magn. Mat. 262, 258 (2003). 

9. F. Mancini, F. P. Mancini, Condens. Matter Phys. 11, 543 
(2008). 

10. K. Maisinger, U. Schollwock, S. Brehmer, H. J. Mikeska, 
S. Yamamoto, Phys. Rev. B 58, R5908 (1998). 



where 



a = e-P v -l = K-l, 
d = e- 2 ' 3V - 1 = K 2 - 1. 



(47) 



We conclude this Appendix by reporting the analytical 
solution obtainable at half filling. Particle hole symmetry 
requires that at half filling the chemical potential must 
take the value /i = U/2 + 2V. For this value of /i, Eqs. 
(|4"5j) become 

Xi = 2GK- 2 {1 -Xi- dX 2 )(l + aXi + a 2 X 2 ) 

+ K-\2 + (d - l)X x - 2dX 2 ]{l + dX x + d 2 X 2 ), 
X 2 = K- 4 [l + dX t - {2d + 1)A 2 ](1 + dXi + d 2 X 2 ) 



2GK- 2 {1 + d)X 2 (l + aX x + a 2 X 2 ), 



(48) 



where G — I" 2 cosh(/3/i). The solution of the above 
equations is: 



X\ = \- 
X 2 



(K + !)[! + 2GK + K 2 - Q] 
AGK{K - 1) 
1 + 2GK + K 2 -Q 
AGK{K - l) 2 ' 



(49) 



